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Abstract - Non-Abelian states are characterized by the existence of zero-mode Majorana fermions 
bound in the quantized vortices, which in the bulk is characterized by an odd Chern number. In 
this Letter, we find that in topological semimetals with a single two-band-touching node all the 
gapped pairing states are non-AbeUan states, regardless of the pairing details. The property is 
further generalized to situations where such two-band-touching are gapped and/or deformed. The 
nonzero Berry phase on the Fermi surface plays an important role in the topological properties of 
the pairing states. 
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O ■ Introduction. — Since the discovery of the quantum 
^~Hall effects [T], efforts devoted to understanding various 
f — .topological states of matter and their phase transitions 

J> greatly enriched the basic concepts of condensed mat- 
'^ter physics [2H3]. One of the significant aspect is that 






X 



two topologically distinct states, which can not be adi- 
_il abatically connected to each other, can have the same 
•symmetry and hence breaks down the Landau-Ginzburg 
L^ paradigm of phase transitions. Moreover, exotic excita- 
jtions obeying non-Abelian statistics have been found in 
,— I genuine and model condensed matter systems [SJIS]. In 
L|6ne known paradigm, the non-Abelian states are closely 
related to superconducting/superfluid pairing states with 
odd Chern number [7], in which there is a topologically 
protected Majorana fermion in each quantized vortex. 
The braiding of Majorana fermions gives rise to the non- 
Abelian statistics [5]. The non-Abelian statistics is pro- 
posed to be exploited to realize the topological quantum 
computation [9HTl]. 

Besides the known non-Abelian states in fractional 
quantum Hall system [5], spin liquids 0, ^He [H] and 
Sr2Ru04 [13], recently there are theoretical proposals of 
realizing non-Abelian states in spin-orbit coupled semicon- 
ductor heterostructures [T5HT5] and in cold-atom systems 
[13]. Especially, the proposal of achieving non-Abelian 
states on the surface of topological insulators proximity 
to s-wave superconductors invoked a lot of studies [15] . 



Semiconductor nanostructures with strong spin-orbit cou- 
pling proximity to superconductors are thought to be an- 
other good candidate to realize the long persued Majorana 
fermions and topological quantum computations [T51IT5] . 
Inspired by those studies in this paper we search for non- 
Abelian states in topological semimetals and related sys- 
tems due to superconducting proximity effect. 

Topological semimetals are systems in which there are 
band-touching nodes with nonzero winding number. In 
particular, here we focus on the situation with single two- 
band-touching node. The winding number is related to 
the Berry phase enclosed by any paths enclosing the node, 
N^ = j^ f^dk- (*(k)|iVk|*(k)) where C is a contour en- 
closing the two-band touching point K in k-space and 
\l/(k) is the wavefunction (single- valued and continuous 
function of k) in the band above (or below) the node. 
Such nodes can be viewed as k-space vortices [THj. Ex- 
amples are Dirac cones and quadratic band touching [see 
Fig. 1], where the winding number is iV„ = ±1 and ±2 
respectively. When the Fermi level is away from the two- 
band touching node, the Fermi surface encloses the node 
and gains a nonzero Berry phase. The Berry phase on 
the Fermi surface plays an essential role as it modifies the 
effective pairing between paired states and hence the topo- 
logical property of the superconductor. 

Remarkably, we find that, in topological semimetals 
with a single k-space vortex, all gapped pairing states 
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are non-Abelian states, regardless of the specific details 
of pairing. We further generalize such property to situ- 
ations where such two-band-touching are gapped and/or 
deformed. The paper is organized as follows: we first dis- 
cuss topological semimetals with a k-spacc vortex carrying 
even winding number. We then study the situation when 
the winding number is odd, which can be thought as gener- 
alizations of single Dirac cone systems. As we aim at very 
general properties of the system, in all the studies, we con- 
sider a general form of pairing, i.e., we assume all the s, p, 
d ... -wave pairing interaction coexist, in contrast that in 
the literature mainly s-wave pairing interaction is consid- 
ered. We then discuss the situations when i) the k-spacc 
vortex is gapped and ii) both gapped and deformed so that 
the two bands go in the same direction. The situation ii) 
for systems with odd winding number can be viewed as 
generalizations of spin-orbit coupled two-dimensional sys- 
tems [IB] . 

Topological semimetals with a single k-space vor- 
tex carrying even v^rinding number. In spinless 
(or spin-polarized) many-fcrmion systems in 2D lattices 
with multiple orbits in each unit cell, there may exist two- 
band-touching nodes (k-space vortices) in the spectrum. 
In systems with a single k-space vortex as well as time- 
reversal and space inversion symmetry, the general form 
of the Hamiltonian is 



ilo(k) = ho{k)ao + /i:,(k)CT^ + h^{k)a. 



(1) 



Here the Pauli matrices act on the Wannier orbits (pseudo- 
spins), (To is the 2x2 identity matrix. Due to time- 
reversal symmetry, the two-band touching node can only 
be at a time-reversal invariant momentum K with — K = 
K as there is only a single node. In the above equa- 
tion, k is measured from K. /i^(— k) = /ii/(k) for 
ly ~ 0,x,z and /ij,(k) = due to the time-reversal 
and space inversion symmetry. The spectrum is ek± = 
/io(k) ± Jhl(k) + hl(k). For semimetals, |/io(k)| < 



ft.^(k) 4- /i^(k). The eigenstates of -ffo(k) are 



l"+(k)) = 2 



l"-(k)) 



-i<t>k 



_e-«0k _ 1 



(2) 



with (/)k = Arg[/iz(k) -I- ihxCk)] being the direction of the 
pseudo-spin polarization in z-x plane. The winding num- 
ber of the k-space vortex can also be written as 



N„ 



1 

2^ 



(3) 





Fig. 1: (Color online) Illustration of the (a) quadratic band 
touching (b) Dirac cone systems with Fermi surface above the 
touching point. Illustration of the direction of the field {hx, hz) 
or {hx, hy) (also represents the pseudo-spin or spin polarization 
direction) at Fermi surface for (c) quadratic band touching and 
(d) Dirac cone system. 



One example of such two-band touching is the quadratic 
band touching in the checkerboard lattice systems [2011^ . 
It is shown [21] [22] that in such systems in the vicinity 
ofK- (^,7r), /io(k) =toe, /i,(k) 



tz[k^ 



kl) 



^ixf^xf^yi hizy^) 



and /iy(k) = with Iq, tx, and tz (l^ol < 



l^xM^zl) being the band parameters. When the Fermi 
level is away from the band-touching point K, there is 
only one band crosses the Fermi level. In this situation, 
the wavefunction (or the pseudo-spin polarization) on the 
Fermi surface winds N^ times, as which encloses the band- 
touching node K. 

The general form of the Bogoliubov-de Gennes (BdG) 
Hamiltonian \s H ^ i X^k *^(k)^k'J'(k), where *(k) = 

(^.t(k), Vi(k), Vr(-k), ^'^(-k)) , and 



Hk = 



ffo(k) - M 

A*(-k) 



-A(k) 
-ffo(-k) - 



A* 



(4) 



Here A(k) = iAo(k)o-,y + A3(k)o-^ + iAy(k)CTo - A:r(k)(T^ 
with Aq and A^^ {v = x, y, z) being the singlet and triplet 
pairing respectively. The above Hamiltonian represents 
a general form of pairing, i.e., includes all the s, p, d ... 
-wave pairing interaction. 

In the weak pairing regime where lAj^l <C |/^|, one 
can ignore pairing between states separated > |/x|. As 
there is only one band crosses the Fermi level and an- 
other is well below, one can project into the subspace 
with only the band crossing the Fermi level. To the lead- 
ing order, the projected BdG Hamiltonian is -ffpBdG 
iEk*p(k)H^*p(k)with*p(k) = (ck±, ■ 



J 



-\l±I 



and 



which is physically transparent: N^ is the winding of the 
direction of the pseudo-spin polarization (or the direction 
of the field h) on the contour C. As /i,y(— k) = /ii/(k), the 
winding number can only be an even integer. 



rii, 



£k± -M 

A:ff(k) 



Aeff(k) 

-£k± + M 



(5) 



Here the -I- and — indices are applied to the /i > and 
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< cases respectively. The effective pairing is 



Aeff(k) = e*' 



iAy - -sgn(/i) ^(A, ± zA,)e±'' 



(6) 

The eigenstates can be obtained by directly diagonalizing 
the above Haniiltonian. The Chern number is given by 
[23] 



Nc 



E- 



rfke^-[VkX (*„|iVk|*„)], (7) 



where ^„ are the wavefunctions of the occupied bands. 
Direct calculation yields [for details, see Appendix] 



A^< 



c 



?ii(/i) 



27r 



|^a.,0A(k) 



(8) 
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where 6'A(k) = Arg[e-"^''Acff(k)]. The Chern number 
is nothing but the winding number of e"*"^'' Aeff (k) at the 
Fermi surface (denoted as 'FS' above). The Chern number 
is thus fully determined by the effective pairing at Fermi 
surface, which is due to the fact that the superconducting 
gap is opened at the Fermi surface [5J. From Eqs. (O and 
([S]), the Berry phase contribution to the Chern number, 
through e^"^'', is clearly seen. 

The error of the eigenstates obtained from the projected 
Hamiltonian is on the order ofOdDj/l/xl). However, this 
induces no error in the calculated Chern number of the su- 
perconductor. This is because one can adiabatically tune 
all the pairing Ai^ (y = 0, x, y, z) to infinitesimally small 
by a universal scaling, which does not close the gap and 
hence keeps the Chern number unchanged. In such a way 
the error can be tuned to infinitesimally small. Therefore 
the above approach faithfully deduces the Chern number 
in the weak pairing regime. 

A crucial observation is that e^"^''Aoff(k) is always an 
odd function of k as A^ (i/ ~ x, y, z) are odd functions 
and e^^'^^ are even functions of k. Hence the Chern num- 
ber Nc is odd for all the gapped pairing states (for nodal 
states the Chern number is ill-defined) . Therefore all the 
gapped pairing states arc non-Abelian states, regardless of 
the specific pairing interactions. 

Topological semimetals with a single k-space vor- 
tex carrying odd winding number. In spinful 2D 
many-fcrmion systems there may exist two-band touching 
points in the spectrum as k-space vortex. The fermion 
doubling theorem states that in 2D lattice systems there 
can only be even number of such vortices [33]. How- 
ever, in the boundaries of two three-dimensional systems 
with different Z2 topology, there can be odd number of 
such k-space vortices. Such as a single Dirac cone on 
the surface of strong topological insulator. The systems 
concerned here have a single such k-space vortex at K 
(with — K = K) and time-reversal symmetry. Rather than 
N.uj = ±1, here we assume the winding number can be any 



odd integer. This is reasonable as the Z2 topology guaran- 
tees odd number of Dirac cones on the surface [1] , which if 
merged by adiabatic tuning of the Hamiltonian can form 
a single k-space vortex carrying any odd integer winding 
number. The general Hamiltonian of such k-space vortex 
is 

Ho{k)^ho{k)ao + h{k)-a, (9) 

where the Pauli matrices now denote true-spin and k 
measured from K. It is required that \ho\ < |h| so 
that the k-space vortex is well-defined. Without chang- 
ing the winding property of the Dirac cone, one can set 
/iz(k) = (This can also be achieved by adiabatic tun- 
ing of the Hamiltonian of the system.). The spectrum is 

e_j_i(; = /io(k)±W/i2 -I- /i2 and the eigenstates are \u±(k)) ~ 

^(e-^'^H t) ± I i)) with Vk = Arg[/i,(k) + ihyik)] is the 
direction of spin polarization. 

The winding number of such k-space vortex can also be 
written as 

^^ = ^i ^^k. (10) 

Again it is transparent that A^^; is the winding number 
of the direction of the spin polarization (or the direction 
of the spin-orbit field). The winding number can only be 
an odd integer as h(— k) = — h(k) due to time-reversal 
symmetry in spinful system. In particular, N^ = ±1 for 
Dirac cone. 

Following the argument in previous section, in the weak 
pairing regime one can study the topological property of 
the system with the projected BdG Hamiltonian [Eq. ([S])]. 
But here 



Aeff(k) 



o^-^k 



1 



sgn(^)Ao-f-^(A,T*A,)e±^'''- . (11) 



The Chern number Nq is the winding number of 
e~*'''''Aoff (k) at the Fermi surface as in Eq. ([5]). It is noted 
that e~*''^''Acff (k) is always an even function of k. There- 
fore the Chern number of each gapped pairing state is 
even. 

As the concerned system is on the boundary of two 
topologically distinct three-dimensional systems, it does 
not have well defined edges [1]. One way to circumvent 
this problem is to circulate the superconducting state with 
a ferromagnetic insulating state with the same ^o(k) as 
well as a magnetization along z-direction Mcr^ (TS] ■ When 
\M\ > 1^1 [26], the quasi-particles can not propagate into 
the ferromagnetic region. On the boundary between the 
superconducting region and the ferromagnetic one, there 
are gaplcss Majorana edge states. The ferromagnetic insu- 
lating state is topologically equivalent to a superconduct- 
ing massive Dirac fermion system with \M\ > |^|. It has 
a Chern number of sgn{M)Nw For instance, there may 
be Uc clockwise moving edge states and Ua anti-clockwise 
moving edge states. According to bulk-edge correspon- 
dence, Uc — Ua = Nc — sgn{M)Nj^ . The difference Uc — n^ 
is fixed by topology and is always odd as Nc is even and 
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N^ is odd. Therefore the total number of edge states 
■AAedge = nc + Ua is definitely odd. 

The above analysis can also be applied to the Majo- 
rana bound states in the core of a quantized vortex, which 
can be viewed as edge states live in the small circular 
edge of the vortex with vacuum at the center j7| . As the 
boundary condition at the center does not affect the ex- 
istence of the zero-energy Majorana bound state, it can 
be tuned that the vacuum at the center is a supercon- 
ducting massive Dirac fermion with \M\ > \^\. Therefore 
there are A/'edgc number of Majorana states in core of a 
quantized vortex. In reality, there are inevitable mixing 
between those states (due to, e.g., disorder) or interac- 
tion between them, which lift the degeneracy. However, 
the particle-hole symmetry guarantees the existence of one 
zero-energy Majorana bound state if ACdgc is odd. Accord- 
ingly, all the gapped pairing states are non-Abelian states 
as A/'edgc is definitely odd. This argument (essentially the 
same as that in Ref. [7]), verifies the existence of the Ma- 
jorana zero modes in the vortex without explicitly solving 
the Schrodinger equations as such property is dictated by 
the bulk topology j7]. Recent theories also present dif- 
ferent proofs of the relation of the number of Majorana 
zero modes to the Chern number (difference) [57] which is 
consistent with the discussion here. 





Fig. 2: (Color online) Illustration of the dispersions in (a) 
gapped (b) gapped and deformed quadratic band touching sys- 
tems as well as (c) gapped and (d) gapped and deformed Dirac 
cone systems. 

Systems with a single gapped/deformed k-space 
vortex. — The general Hamiltonian is 

i7o(k) = /io(k)<To -I- h^{k)a^ + h,{k)a, + hy{k)ay. (12) 



The spectrum is £k± = ho{k) ± Jhl{k) + hl{k) + hl{k). 
When k-space vortex carrying even winding number, the 
eigenstates are 

1 



l«+(k)> = 



cos^e 
icos^e" 

I sm %e" 



sin^ 
i sin ^ 



\/2 V - sin ^ 



I cos 

- COS 



2 
2 



(13) 



where ?7k = Arg[/ij^ +iy^hl + hi] and 0k = Arg[/i^ -I- ihx]. 
At fe = 0, /ij: = /i^ = whereas hy ^ 0. At large fc, \hy\ ^ 
tt/2. Wc restrict the discussion to the 



^/h^^+h^ and r^k 

situations where hy{—k) = hy{k). For the k-space vortex 

carrying odd winding number, the eigenstates are 



l"+(k)> = 
l"-(k)> 



cos ^e-*'/'" 

• Ck 
sm ^ 

sin ^e-''^" 



(14) 



ih,, 



where ^k — Arg[/i2 + i^/hl + h^] and ■i/'k ~ Arg[/i2: 

Here at k = 0, h^ ~ hy ~ and h^ ^ 0, whereas at 

large k, \hz\ <^ w/i^ -|- hy and Ck — > 7r/2. We also assume 

h,{-k) = h,{k). 

Let's first consider the case where the single k-space 
vortex carrying even winding number is gapped. When 

\hy{0)\ < \fi\ and \ho\ < J hi + h"^ + h^ [see Fig. 2a] one 

can similarly reduce the problem to that described by the 
projected BdG Hamiltonian in Eq. ^ for each band in 
the weak pairing regime. The effective pairings are 



Aelf(k) 



eJ'l"'(±iAySmr]k 

— Ax [cos 0k ± i sin 0k cos J^k] 

— A2[±icos0kCOS77k — sin0k] f, (15) 

where -I- and — are for the higher and lower bands respec- 
tively. Direct calculation yields that the Chern number is 
still given by Eq. ([8]), which is because a gap opened below 
or above the Fermi level can be closed by adiabatic tun- 
ing. A nontrivial situation is when the k-space vortex is 

both gapped and deformed so that ho > \ h^. + h"^ + hi . 
In this situation the two bands go in the same direction 
at large k [see Fig. 2b]. In this case, when \^\ < \hy{0)\, 
following the derivation in Appendix, one can show that 
the Chern number is 



Nr 



•IlT 






^A(k) 



Wk^A 



sgn[hy{0)]N,, 



(16) 
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where 9^ = Arg[e~''/"'A^fj(k)]. It is seen that as N^ is 
even, the Chern number is again odd for all the gapped 
states. Therefore in this situation all the gapped pair- 
ing states are non-Abelian states regardless of the specific 
pairing details. When the Fermi level is so high such that 
both bands cross it, the total Chern number consists of 
the contributions from each band, which is 



Nc 



E 

± 



2tt 



■deJ 



A(k) 



(17) 
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Hence the total Chern number becomes even (Abelian 
topological superconductor) , when the Fermi level crosses 
two bands. 
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Now we turn to the situation where the single k-space 
vortex carrying odd winding number is gapped. The effec- 
tive pairing in the two bands are 



Aelf(k) = 



- e'^f"^ 



+ sin^(^)(A^T«Aj,)e'" 



cos2(^)(A,±iAj,)e' 



2 



^±ii>k 



T Ao sin(Ck) 
(18) 



Again the Chern number is the same as that at hz — 

when \fi\ > \hz{0)\ and |/io(k)| < J hi + h^ + hj [see 

Fig. 2c] as the gap is opened below or above the Fermi 
level. When the k-space vortex is gapped and deformed, 



so that /io(k) > \/h1 + h"^ + hi where the two bands go in 
the same direction at large k [see Fig. 2d], at |/x| < |/ix(0)|. 



one can show that 



Nr 



27r 



deJ-^0^) 



sgn[;i,(0)]iV^ 



(19) 
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with e% = Arg[e"*'''''A^ff(k)]. It is seen that e-*'^''Acff(k) 
is always an even function of k as e^^"^^ , A^, and Aj^ are 
odd functions while Ck and Aq are even functions of k. 
Thus the Chern number is odd for all the gapped pairing 
states. When the Fermi level is higher so that both bands 
cross it, the total Chern number is given by Eq. (J17p . which 
is always even. 

Relevant physical systems. — Besides the systems 
already discussed in the literature, such as systems with 
single Dirac cone and semiconductor quantum wells with 
Rashba spin-orbit coupling, there arc many unexplored 
candidate systems which fit into the above discussions. 
Below we list some examples of physical systems that re- 
alize the models discussed above: 



Conclusion and discussions. — In this paper 
we discussed the topological properties of two classes of 
semimetal systems under superconducting proximity ef- 
fect in general situations. Those two classes are semimet- 
als with a single two-band-touching node which carries 
i) even ii) odd winding numbers respectively. Remark- 
ably, we found that in those systems all the gapped pair- 
ing states are non-Abelian states regardless of the pairing 
details. This property is further generalized to systems 
when such two-band touching is gapped and/or deformed. 
Those findings supply useful information and a route for 
the search for non-Abelian states. 

Although some of special systems, such as single Dirac 
cone system }15j and two-dimensional electron system with 
Rashba spin-orbit coupling [16j . have been discussed in 
the literature (results here are consistent with those), the 
results obtained in this paper reveals a more general pic- 
ture of the topological property in a general topological 
semimctals and related systems with general pairing inter- 
actions. The theory presented above can be directly gen- 
eralized to systems with multiple two-band touching nodes 
when the pairing between different nodes can be neglected. 
It is noted from Eq. ^ that the Chern number changes 
sign when the chemical potential goes through the band- 
touching node, which indicates that there is a topological 
phase transition in the strong pairing regime. From spon- 
taneous symmetry breaking under attractive interaction, 
the gapped pairing states usually gain more condensation 
energy compared to the nodal ones and thus become en- 
ergetically favored |13]. Hence the above special property 
also implies that the non-Abelian states may be energeti- 
cally favored to be the ground states in the superconduct- 
ing/superfluid phase [22j . 



• Semiconductor nanostructurcs with Zeeman (ex- 
change) sphtting and arbitrary spin-orbit coupUng. 
Given that the winding number at Fermi surface 
is nonzero, two-dimensional systems arbitrary spin- 
orbit coupling and Zeeman splitting can achieve non- 
Abelian states when there is only one band crossing 
the Fermi level. This is a direct generalization of the 
studies in the literature. 

• Thin films of topological Weyl semimetals. In 
Rcf. HH], it is found that in the thin film of topo- 
logical Weyl semimetal HgCr2Se4 the Chern number 
depends on the thickness of the film. At the critical 
thickness, there will be a single two-band touching 
with even winding number N^j — ±2. Around the 
critical thickness, there are gapped quadratic band 
touching which can also give rise to triplet pairing 
and non-Abelian states. 

• Optical lattices with a single quadratic band touch- 
ing. Such as, checkerboard lattices and kagome lat- 
tices noilll]. 
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Appendix: Details of the derivation of the Chern 
number. — Consider, e.g., systems with a single k- 
space vortex carrying even winding number with fi > 0. 
There are two occupied bands after diagonalizing the BdG 
Hamiltonian: one comes from the band crossing the Fermi 
level, the other from the band below the chemical poten- 
tial. In the weak pairing regime where |Ai,| <C |/i|, one can 
ignore pairing between states separated > |^|. By doing 
so, one can obtain the approximate wavefunctions of the 
former and latter occupied bands, which are: 



/si 



sm % cos ^e'"^ 



*o 



^ g»0k/2 



\ 



V -^ 
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*„(k) = e-*'^''/2 



/ sin^ \ 

-cos% 



V / 



(20) 



respectively. Here ^k = Arg[ek+ ^ /^ + «|Acff(k)|] and 
0A = Arg[e~*'^''Aeff(k)]. An important property is that 
the Chern number Nc does not change without closing 
the gap. One can then simplify the calculation of Nq 
by adiabatically tuning the system. It is noted that the 
gap is determined by |Aeff(k)| at the Fermi surface. One 
can then adiabatically tune the system so that |Aeff(k)| 
is nonzero only in the vicinity of the Fermi surface pi] . 
The angular dependence |Aeff(A:,0k)| (here kx = kcosOk 
and ky = k sin 9^) at each energy contour can also be 
adiabatically tuned to be identical to that on the Fermi 
surface. The Chern number is the integration of the Berry- 
curvature over the first Brillouin zone, 



^c=J2^ I dke. 



[Vk X (*„KVk|*r 



(21) 



One can divide the contribution of the integration into two 
parts: one from integration over small k (with a cut-off A), 
another from integration over large k region. Since there 
is no band-gap closing in the large k region, the Chern 
number is determined in the small k region, especially, in 
the vicinity of the Fermi surface where the pairing gap 
evolves, as shown in Ref. [24]. Inserting Ea. ((20|) . one can 
show that the Chern number due to the v band is zero. Nc 
is then solely determined by the o band. Direct calculation 
yields 



Nc 



1 

'2^ 
1 



dOu / dkdk 
Jo 

2tt 



sin2|(a,,0A) 



= 7^ / dOkdgjA 

ZTT 



(22) 



FS 



which results in Eq. ([8|). Other equations for the Chern 
number in the main texts can be derived similarly. 

REFERENCES 

[1] K. V. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 

45, 494 (1980); D. C. Tsui, H. L. Stormer, and A. C. 

Gossard, Phys. Rev. Lett. 48, 1559 (1982). 
[2] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. 

den Nijs, Phys. Rev. Lett. 49, 405 (1982). 
[3] X. G. Wen, Quantum Field Theory of Many -Body Systems 

(Oxford University Press, Oxford, 2004). 
[4] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010); X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 

1057 (2011). 
[5] G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 
[6] A. Kitaev, Ann. Phys. 321, 2 (2006). 
[7] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
[8] D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 
[9] C. Nayak et al. Rev. Mod. Phys. 80, 1083 (2008). 



[10] 

[11] 
[12] 
[13] 
[14] 

[15] 
[16] 



[17; 

[18 

[19 
[20 

[21 

[22' 
[23 

[24' 

[25 
[26 

[27; 

[28; 

[29; 



M. Freedman et al.. Bull. Am. Math. Soc. 40, 31 (2002); 
A. Y. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 
S. Tewari, S. Das Sarma, C. Nayak, C. Zhang, and P. 
Zoller, Phys. Rev. Lett. 98, 010506 (2007). 
G. E. Volovik, The Universe m a Helium Droplet, (Claren- 
don, Oxford, 2003); arXiv:1111.4627 

G. M. Luke, et al, Nature 394, 558 (1998); A. P. Macken- 
zie and Y. Maeno, Rev. Mod. Phys. 75, 657 (2003). 

C. Zhang, S. Tewari, R. M. Lutchyn, and S. Das Sarma, 
Phys. Rev. Lett. 101, 160401 (2008); M. Sato, Y. Taka- 
hashi, and S. Fujimoto, ibid. 103, 020401 (2009); S.-L. 
Zhu et al, Phys. Rev. Lett. 106, 100404 (2011). 

M. Sato, Phys. Lett. B 575, 126 (2003); L. Fu and C. 
L. Kane, Phys. Rev. Lett. 100, 096407 (2008); A. R. 
Akhmerov, J. Nilsson, and C. W. J. Beenakker, Phys. 
Rev. Lett. 102, 216404 (2009); C.-K. Lu and I. F. Her- 
but, Phys. Rev. B 82, 144505 (2010). 
M. Sato and S. Fujimoto, Phys. Rev. B 79, 094504 (2009); 
J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 
Phys. Rev. Lett. 104, 040502 (2010); J. Alicea, Phys. Rev. 
B 81, 125318 (2010); R. M. Lutchyn, J. D. Sau, and S. 
Das Sarma, Phys. Rev. Lett. 105, 077001 (2010); Y. Oreg, 
G. Refael, and F. von Oppen, ibid. 105, 177002 (2010); 
L. Santos, T. Neupert, C. Chamon, and C. Mudry, Phys. 
Rev. B 81, 184502 (2010); M. Sato, Y. Takahashi, S. Fu- 
jimoto, ibid. 82, 134521 (2010); L. Santos, Y. Nishida, C. 
Chamon, and C. Mudry, Phys. Rev. B 83, 104522 (2011); 
L. Mao, J. Shi, Q. Niu, and C. Zhang, Phys. Rev. Lett. 
106, 157003 (2011). 

X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 
82, 184516 (2010). 

J. Alicea, Y. Oreg, G. Refael, F. von Oppen and M. P. A. 
Fisher, Nat. Phys. 7, 412 (2011). 
G. E. Volovik, Lect. Notes Phys. 718, 31 (2007). 
K. Sun, H. Yao, E. Fradkin, and S. A. Kivelson, Phys. 
Rev. Lett. 103, 046811 (2009). 

K. Sun, W. V. Liu, A. Hemmerich, and S. Das Sarma, 
Nat. Phys. 8, 67 (2012). 
J.-H. Jiang, arXiv:1110.2043 

D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 
1959 (2010). 

X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 

81, 134508 (2010). 

H. Nielssen and N. Ninomiya, Phys. Lett. 130B, 389 

(1983). 

More precisely |Af| > yTi^ + Af as the s-wave pairing 

modifies the energy spectruin of the Bogoliubov quasi- 

particle at k = [TS] . 

L. Santos, Y. Nishida, C. Chamon, and C. Mudry, Phys. 

Rev. B 83, 104522 (2011); R. Roy, Phys. Rev. Lett. 105, 

186401 (2010). 

G. Xu, H. Weng, Z. Wang, X. Dai, and Z. Fang, Phys. 

Rev. Lett. 107, 186806 (2011). 

P. W. Anderson and P. Morel, Phys. Rev. 123, 1911 

(1961); M. Cheng, K. Sun, V. Galitski, and S. Das Sarma, 

Phys. Rev. B 81, 024504 (2010). 



p-6 



